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Abstract. - The properties of the Kohn-Sham (KS) exchange potential for open systems in 
thermodynamical equilibrium, where the number of particles is non-conserved, are analyzed 
with the Optimized Effective Potential (OEP) method of Density Functional Theory (DFT) 
at zero temperature. The quasi two-dimensional electron gas (2DEG) is used as an illustrative 
example. The main findings are that the KS exchange potential builds a significant barrier-like 
structure under slight population of the second subband, and that both the asymptotic value 
of the KS exchange potential and the inter-subband energy jump discontinuously at the one- 
subband (IS) — > two-subband (2S) transition. The results obtained in this system offer new 
insights on open problems of semiconductors, such as the band-gap underestimation and the 
band-gap renormalization by photo-excited carriers. 



Density Functional Theory (DFT) has become a most used computational tool in the study 
of strongly inhomogeneous interacting systems such as atoms, molecules, clusters, and solids. 
[1] A crucial step for any DFT calculation is the approximation used to evaluate the exchange- 
correlation (xc) contribution E xc = E x + E c to the total energy. Traditionally E xc has 
been evaluated by using density-dependent ic-functionals, either of local (LDA) or semi-local 
(GGA, meta-GGA,...) character. [2] In the last few years, however, considerable interest has 
arisen around orbital-dependent xc functionals, which are implicit functionals of the density 
and whose implementation in the KS scheme is known as the OEP method. [3] If the Hartree- 
Fock expression for the exchange energy functional is used, and E c is neglected, the OEP 
method is equivalent to the exact x-only implementation of KS theory. [4] Several advantages 
are associated to the use of exact x-only DFT calculations for (closed) systems with a fixed 
number of particles: cancellation of the spurious Hartree self-interaction energy, [5] correct 
high-density limit, [2] great improvement in the KS eigenvalue spectrum, [6] semiconductor 
band structure and excitations, [7] and nonlinear optical properties. [8] It is the aim of this 
work to report on further novel features of the OEP scheme for open configurations where 
the system under study is in thermodynamical equilibrium with a reservoir and can exchange 
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particles with it, using a quasi- two-dimensional electron gas at zero temperature to illustrate 
our results. We compare usual approximations against our exact (OEP) x-only results. 

The present theory can be applied to any system with translational symmetry in a plane. 
For instance, a semiconductor modulation-doped quantum well (QW) grown epitaxially as 
shown in the inset of Fig. ^ In these systems, it is possible to confine an electron gas 
by changing the semiconductor in the growth direction z. If the larger gap semiconductor is 
doped with donors, it provides electrons to the trap formed by the smaller gap semiconductor, 
until the charge-transfer field equilibrates the common chemical potential /i in the QW and 
a doped region acting as particle reservoir. The charge transferred from the reservoir to 
the QW can be tuned by an external electric field. Assuming translational symmetry of 
the 2DEG in the x — y plane (area A), and proposing accordingly a solution of the type 
<^ker(p, z) — exp(tk-p)£f (z)/VA, the ground-state electron density can be obtained by solving 
a set of effective one-dimensional Kohn-Sham equations of the form 

where effective atomic units have been used. £f (z) is the eigenfunction for electrons in subband 
i (= 0, 1, ...), spin a (=|, |), and eigenvalue sf. The local (multiplicative) Kohn-Sham potential 
V£ s (z) is the sum of several terms: V^ s (z) — V ext {z) + Vh(z) + V£ c (z). V ext {z) is given 
by the sum of the epitaxial potential plus an external electric field. Vh{z) is the Hartree 
potential. Within DFT, V£ c (z) — 8E XC / 8n a (z). In an open system, E xc = E xc [{e% }, {£f (2)}] 
is a functional of the set of ef's and £f (z)'s. The zero-temperature 3D electron density is 

n a {z) = £ e « <M (fcj?) 2 l£ T ( z )| 2 / 47r ; with k F = V 2 (P - £ D- After some len S th y but standard 
manipulations of the OEP scheme, [3,9] the calculation of V xc {z) for an open system can be 
summarized in the following set of equations: 

occ. 

i 

S ia (z) = (k$) 2 W(znf(z) - CZ \%{z)\ 2 + c.c, 

The sum in Eq.© runs over the occupied subbands, C% = V l xc +{2ir/A)dE xc /de° , AV™ (z) = 
V° c (z) - u%(z), u™ c {z) = 4:TT/A(k l £) 2 tf (z)*5E xc /6£? (z), and mean values are defined as 

O = J dztf(z)*0(z)g(z). Eqs. (2® dete rmine the local V xc (z) corresponding to an or- 
bital and eigenvalue-dependent approximation for E xc . Eqs. Q an< i © have to be solved 
self-consistently. 

Some comments are necessary here: a) Our treatment corresponds to an open system, and 
this is the origin of the term V^. in C l xc on the r.h.s. of Eq.(J3J. This term gives a finite (7^ 0) 
contribution under the replacement V xc (z) — > V xc (z) + a a , with a a an arbitrary constant. As 
a consequence Eq.JSJ is not invariant if V xc (z) is shifted by a constant, and V xc (z) is strictly 
determined. For closed systems, this term is absent and the corresponding OEP expression 
gives V xc {z) up to an additive constant. [3] b) Using Eq.Q it can be found that the ipi(z) 
(denoted as "shifts") satisfy the orthogonality constraint J dz^f(z)*ipf(z) — 0. Application 
of h^siz) to the 4>?(z) in Eq. (@J yields an inhomogeneous differential equation, which can 
be considered as an alternative definition of the shifts. [10,11] c) Integrating Eq.(j2J and using 
the orthogonality constraint, we obtain the important property YIT = 0- 



2dz 



(1) 



(2) 

(3) 
(4) 
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Using Eqs.Q and P)l. and dropping spin indices by assuming a paramagnetic situation, 
V xc (z) could be formally expressed for real £i(z)'s (as is our case) as V xc (z) — V xc i(z) + 
V xc2 {z) + V xcS (z), where 



V Az) = f [kUi(Z)] '' 

v xcl {z) - 2^ 2lTn{z) 



u xc {z) + AV X( 



V xc2 (z) = Y,( k F?Si{z)/[2im{z)}, 

i 

Vxc3{z) = 2Mi) 



(5) 
(6) 
(7) 



with primes denoting derivation with respect to z. It can be shown that if the shifts in V xc (z) 
are forced to be zero, the only term left is V xc i(z). This truncated expression for V xc (z) is 
identical to the one obtained in Ref. [12] for closed systems (KLI approximation). Accordingly, 
we identify here that V^ LI (z) = V xc \{z). In the IS regime V xc2 {z) = V xc §{z) = 0, [13] 
then V^c S \z) = V* LI (z). A KLI calculation in the many-subband situation, requires in our 
open case to find a self-consistent solution of Eq. with V xc (z) — V xc i(z), constrained 
by the exact condition ^2° cc ' C xc = 0. This procedure univocally determines V xc i(z). It 
is important to realize however, that what we call KLI approximation is not the same as 
previously considered [12]. While the expressions for V^ LI are identical in both cases, the 
boundary conditions are different. In our case, we use *}2° cc C xc = 0. In previous studies, 
the physical boundary condition V x ° = u^f was imposed in order to satisfy the correct 
asymptotic behavior (m denoting the index of the highest-occupied subband). All results 
given until this point include both exchange and correlation. In the s-only version of the 
OEP method, V xc -> V x . [4,6,11,12] 

We plot in Fig. Q]the x-only OEP for several fillings. It is noticeable the large difference 
for V x (z) as calculated for A = [i — e-y — * 0~ (A) or A — * + (B), as well as the "hump" 
that V x (z) develops for incipient population of the second subband, and the fast decay of this 
"hump" after macroscopic occupation of the band (C). 

The x-only eigenvalue spectra, in the neighborhood of the IS* — > 2S transition, is shown 
in Fig. El The top panel corresponds to the inter-subband energy e i = e% — eo calculated 
in three different ways: OEP, KLI and LDA. The lower panel displays the asymptotic values 
of the x-only potential. We notice an abrupt increase of £q^ p when the system passes from 
the IS to the 2S regime. This should be contrasted with the abrupt decrease of e^ 1 , and 
the continuity of £qi K - Also noticeable is the abrupt jump of V x (z — > oo) in OEP and KLI at 
A = 0. 

In what follows we will provide a brief explanation for each of the striking results shown in 
previous figures. We start with the asymptotic behavior of V x (z), for the many-subband case 
m > 1. Following the analysis of Ref. [14], and based on the fact that — * oo) — > e~@ iZ 
(disregarding factors involving powers of z) with /3j > and /3j < f3j if Si > Sj, one can derive 
from the asymptotic limit of the differential equation for the shifts that ipi<m(z oo) — » 
e -0 m z_ From the asymptotic analysis of Eq .(01, we obtain ip m (z — > oo) — > [^"/(fc™) 2 ] e~ l3mZ . 
Consequently, all the shifts corresponding to occupied subbands decay exponentially at the 
same rate. Inserting this result in the asymptotic expression for V x (z), which amounts to 
restrict the sum over the occupied subbands to the last one (i = m), it can be checked that 
V X 2(z) and V x ^{z) vanish to leading order, while the contribution from V x \{z) remains finite, 
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Fig. 1 - Thin lines: exchange-only OEP for A = —0.1(A), 0.01(B), 0.1, 0.3 and 10 meV (C). Thick 
line: V x 2(2) in the limit kp — > + (see text below). The arrows on the right denote the asymptotic 
value of the exchange potential for selected subband fillings. Inset: schematic view of our model for 
the modulation-doped QW. The thick vertical stripe on the left represents the ionized donor impurities 
region. The thin vertical line on the right represents a distant metallic plane which induces a charge- 
transfer field along z included in V ex t(z). If the metallic plate is positively (negatively) charged more 
(less) electrons are transferred towards the well. 




Fig. 2 - Exchange-only energy spectra by different methods. Top panel: inter-subband energy spacing 
£01 and fj, — so- Lower panel: asymptotic value of the exchange potential V x (z — > 00) = AV X (see Eq. 
JjjJ). The KLI approximation coincides with the exact result (OEP) in the IS regime. Open circles: 
same as above, but including correlation a la LDA. 
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14(2 - CX>) - V X1 {Z - CX>) -> U™(Z - 00) + ATC, (8) 

where u x n (z — > oo) — > —1/z. [15] Eq.(JHJ coincides with the asymptotic result found in standard 
applications of the OEP method to closed systems. [3] In such discussed above, 

V x (z) is determined up to an additive constant, which is fixed arbitrarily choosing AV X 
0. In our open configuration, however, the constant is automatically generated by the self- 
consistent solution of Eqs.([p)-I0). Remembering that AV X — V x — m™, we have no reason to 
expect a smooth change in the value of the asymptotic constant at the IS* — > 25 transition. 
This explains the abrupt jump at A = in Fig. The KLI approximation preserves this 
jump in AV X , but fails quantitatively, giving a jump a 23 % larger than the exact one. 
V x BA (z — > oo) = 0, independently of subband filling. It is worth at this point to address 
the important difference with the related results of Refs. [6, 11]. These works are concerned 
with the asymptotic behavior of the x-only OEP V x (r) as applied to closed systems such 
as isolated molecules and metallic clusters. Their main finding is that V x (r) can approach 
different asymptotic values depending on the spatial direction with which the asymptotic 
region is reached. In contrast, the solution of this open system shows abrupt changes in the 
asymptotic values of the ir-only OEP resulting from changes in the subband filling, at a fixed 
spatial direction (z). 

Let us concentrate now on the the "hump" effect observed in Fig. ^ This effect is due to 
the contribution V X ci{z) in V xc (z). To see this, let us analyze V XC 2(z) in the limit kp — > + . In 
this almost depleted 25 regime, we can approximate V XC 2(z) — C xc — (k F ) 2 i/)o(z)/£o(z). Using 
S (z) + S 1 (z) = 0, and C° c + C xc = 0, we find that Vo(^) - C° xc [^{zf - ^(z) 2 ] /£o(z)(k F ) 2 . 
Replacing ifio(z) in the latter expression for V xe 2(z), we obtain V XC 2(z) ~ C xc [£i(z) / (,o( z )] 2 ■ 
We plot this approximation for V X 2(z) in Fig. ^ (thick line), showing that it is the "limiting" 
value of the full V x (z) if A — > + in the region where the ratio £i{z)/£o(z) is large, excluding 
the asymptotic region, which is dominated by V x \{z). Since V xc i (and V xc z) are neglected in 
the KLI approximation (as defined above), the "hump" effect is also lost. From this point of 
view, our KLI is not a good approximation for slight occupancies of excited subbands. 

To further clarify the role played by the different contributions to the x-only OEP, we 
display them in Fig. [3] It is clear from this figure that the "hump" effect is due to V x z{z), 
while V x i(z) dominates in the high-density region and gives also the main contribution to 
V x (z) in the asymptotic region, according to the result of Eq.(JH|. The contribution from 
V x s{z) has been found to be small in all the analyzed situations. 

The building of this significant barrier explains, in turn, the abrupt increase of £o\ EP under 
second subband occupation: as the barrier lies about the same zone where the weight of 
£i(z) is concentrated, it affects mainly the electrons in the second subband, resulting in an 
increase of the inter-subband spacing £o EP . However this increase is the result of two opposite 
contributions: the barrier building, which tends to increase £oi EP , and the abrupt jump in 
the value of V x (z — > oo), which tends to decrease £oi EP - [16] The fact that the sum of the 
two contributions result in a net increase of Eqi > implies that the repulsive barrier effect 
overcomes the decrease of V x (z — + oo). This latter effect is the only one included in £qi j that 
leads to an abrupt decrease of e^ 1 at the 15 — ► 25 transition. Finally £qP a , lacking from the 
"hump" and asymptotic value effects, displays a smooth transition at A = 0. Several works 
in the past have addressed the (surprising) issue of the discontinuity of the exact V xc (r) as 
a function of the electron number. [12, 17-19] Such a V xc , for a filled-subshell system with 
N electrons, must essentially jump by a constant for finite r as the number of electrons is 
allowed to vary from N — 5 to N + S where S is a positive infinitesimal. This non-analytic 
behavior of V xc leads in turn to discontinuities on the iV-dependence of the eigenvalues, which 
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Fig. 3 - OEP and LDA s-only potentials for A = 0.3 meV. ^o(z) and £i(z) are the ground and first- 
excited subband wave-functions, and n(z) is proportional to the density. The asymptotic expression 
of V x (z) ~ tt™(z) + AV X (Eq.(8)) is only plotted in the barrier region. V x (full line) and V x including 
correlation effects a la LDA (not shown) are indistinguishable on the scale of the figure. 



are reminiscent of the results presented in Fig. 2. In our case, V x can also in principle change 
discontinuously by a constant at the IS — > 25* transition. Instead, our calculation predicts 
that the system shows another type of discontinuous behavior, replacing V x (z) (curve A in 
Fig. 1) by another function V x {z) (curve B in Fig. 1). The resulting V x (z) induces a charge 
transfer from the minority subband to the other optimizing the exchange energy by maximizing 
overlaps. A constant shift in the exchange potential does not favor such transfer. We have 
assumed that exchange is the dominant contribution although correlations could contribute 
subband discontinuities [19] . In order to estimate the importance of correlation on the x-only 
results, and considering the good agreement between V x and V r x LDA in the high density region 
inside the quantum well displayed in Fig. 3, we show in Fig. 2 with open circles the effect of 
replacing V; OEP (z) by V° EP (z) + V; LDA (z), both on eg EP and V° EP (z -> oo). This amounts 
to re-calculate the self-consistent solutions of Eqs.QJ and (J2J for all A, with V x OEP (z) obtained 
according to the OEP procedure, while V C LDA = 5E EPlA / 5n(z). [9] The correction on £qi EP is 
extremely small (~ 0.1 meV for all values of A), and could be understood as resulting from a 
shift of the x-only value of eoi under a small (perturbative) correlation potential. The impact 
of LDA correlation on the discontinuity at A = in £q EP and V x OEP (z — ► oo) is negligible. 

Let us discuss the results for this system, in the more general context of the semiconductor 
gap problem. [19] Note that if the gap is defined as the energy difference of exchanging an 
electron or hole between a given system and a reservoir [2, 3] then, i) it is necessary to 
consider an open system, and ii) the gap will be affected by a sudden change of V x (r) under 
an infinitesimal occupation of the conduction band as observed in this work. This study 
also suggest that going beyond KLI might be required to describe exchange contributions 
to the gap. The strong decrease of the barrier-like component on V x (z) in Fig. 1 for A ~ 
0.01 meV — > 0.3 meV implies a significant reduction of the exchange contribution to £oi EP . 
This reduction of £oi EP resembles the one observed in the gap of semiconductors under intense 
photo-excitation (band-gap renormalization). The £° EP results of Fig. 2 for A > suggest an 
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exchange contribution to the band-gap renormalization of the same order of magnitude than 
the A = gap discontinuity itself. This should be contrasted with the smaller reduction rate 
with increasing density displayed by e^ 1 and e}^ A . 

In summary, we have shown that the x-only OEP has a number of interesting properties 
when applied to an open 2DEG: i) under depopulation of the highest-occupied subband, V x (z) 
develops a barrier-like structure in the spatial region where the density starts to be dominated 
by electrons in that subband; ii) the size of the barrier increases dramatically when the highest 
subband occupation is tuned to complete depletion; iii) V x (z — > oo) strongly depends on 
subband filling, showing an abrupt jump at the IS* —> 2S transition; and iv) as a consequence 
of the combined effect of ii) and iii), the inter-subband energy is discontinuous at the IS — > 
2S transition. Correlation effects, estimated a la LDA, are found to be very small for this 
system, giving hope for its rigorous and systematic consideration following the lines of the KS 
perturbative approach. [20] The results for the present system offer insights into fundamental 
problems related to the band-gap in semiconductors. 
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